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Abstract. The standard approach to determine the parameters of a resonance is based on the study of the volume dependence 
of the energy spectrum. In this work we study a non-linear sigma model coupled to a scalar field in which a resonance emerges. 
Using an analysis method introduced recently, based on the concept of probability distribution, it is possible to determine the 
mass and the width of the resonance. 
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INTRODUCTION 

Lattice simulations provide a non-perturbative framework for calculating many low energy properties of QCD. Of 
these studies, the mass spectrum of the lightest hadrons is amongst the most important. In a lattice calculation the 
(stable) hadron spectrum is determined by numerical measurement of the energy eigenvalues of the QCD Hamiltonian. 
Nonetheless, most hadronic states are resonances which cannot be identified with a single energy level of the 
Hamiltonian: this complicates their mass and the width determination. A standard approach, introduced by Luscher [1, 
2, 3] (see also [4, 5] and for a generalization in moving frames [6]), to extract resonance parameters, starts by studying 
the volume dependence of the energy spectrum of the Hamiltonian when the system is confined in a finite box. There 
is a connection between the two-particle energy spectrum in finite volume and the elastic scattering phase shift in the 
infinite volume. 

It was noted that where the phase shift in infinite volume quickly passes through n/2 an abrupt rearrangement of 
the finite volume energy levels occurs. This peculiar behaviour, known as an avoided level crossing (ALC), is the 
signature of a resonance in the lattice data. Therefore, it is possible to extract the phase shift from lattice simulations 
and, consequently, to determine the position (the mass value) and the width of the resonances. The feasibility of this 
method has been demonstrated in four dimensional Ising model [7] and for models in two dimensions: in particular 
in the 0(3) non-linear a— model [8] and in two coupled Ising models [9]. Subsequently, it was applied in three 
dimensions in QED to study the meson-meson scattering phase shift [10] and in four dimensional 0(A) <j) 4 theory [11]. 
Recently the method has been used to extract the p resonance from lattice QCD [12, 13] and to calculate the S-wave 
pion-pion scattering length in the isospin 1 = 2 channel and the P-wave pion-pion scattering phase in the isospin 1=1 
channel [14]. 

During the last years many steps have been taken to address the problem of resonances on lattice; however, the 
approach described before is limited. Many of the resonances in QCD are characterized by a large width and so the 
rearrangement of the finite-volume energy levels does not provide a clear signature of the presence of a resonance: the 
ALC is washed out. A possible solution to this problem comes from a new method to analyze the finite volume energy 
spectrum [15]. The idea behind this work is to test this method in the case of a simple field model which is an effective 
low-energy description of QCD. Moreover, we can take advantage of the results in this model to quantitatively predict 
the precision of the results in the determination of parameters of resonances in QCD; we plan to deal with the glueball 
resonance. 



THEORETICAL BACKGROUND 

In lattice field theory we determine the mass of a stable particle from the exponential decay of suitable correlation 
functions in Euclidean space. If G(x,y) is the two-point Green function of a generic field 0, we can consider the 



partial Fourier transform 



C(t,p) = J -0^e-^G(x,t;d,O) = J ^e icot G(p,co) , (1) 

where G(p, (o) is the Fourier transform of G{x 1 y). One then expects that at large times t 

C{t,p = 0)°ce- m O t , (2) 

where is the lightest mass which appears in the spectrum of the theory. On the other hand, if we consider an 
unstable particle which can decay in two particles % (0 — > 2%) when m§ > 2m n , at large times t we have 

C(f,7? = 0)°c e- 2m "' , (3) 

and so we have no information about the mass of m^. 

Theoretically, we can characterize a resonance studing the S-matrix theory; let us consider the element Si(E) (I 
is the angular momentum and E is the energy, considered as a complex variable) then we can define a resonance as 
the pole E of S[(E), in the complex energy plane, when Im^) < and Re(£) ^ 0. In fact, if Im(£) is small it is 
possible to show that the presence of the pole appears as a peak in the total cross section, for the corresponding /-wave 
contribution: 

^ (£) ^ (£ -M)CrV4 (4) 

where M °= Re^) and Y °< lm(E). As a consequence, a resonance is characterized not only by a mass M but also by 
a width F; therefore, we cannot identify a resonance with an isolated energy level. This has to be compared with the 
fact that the energy spectrum of the Euclidean theory is always real, i.e. we can only study isolated levels. It is clear 
that it is not possible to get the parameters of a resonance using the standard tools of the lattice approach. 

As explained in the introduction, if a resonance is present in a theory, it manifests itself as a rearrangement of the 
energy levels related to a mixing between the unstable particle and the two stable particles n; taking advantage of 
this effect, i.e. using the Liischer approach, it is possible to obtain the resonance parameters. 



Interacting particles in a box 

Let us consider a system of two identical non-interacting bosons of mass m n in a box of volume V = nf=i Li- The 
two-particle states are characterized by a total energy E, in the rest frame of the box, and by a relative momentum p 
related by: 

E = 2yjml + \ff. (5) 

The quantized momentum pi are given by /?,• = ^jH;, where n, G Z. Note that Eq. (5) is only an approximation on the 
lattice because it does not take in account space-time discretization artefacts. 

We can study the spectrum of this system using Eq. (5); in particular, we can determine it for different values of 
Li and for different values of n (n — (n\,ti2,n-$)). In a cubic box, if n 2 = £? =1 nf is fixed, we have degenerate energy 
levels for different values of n,-. In Figure 1 (Left) we show the five lowest levels as function of L (cubic box). As 
proposed in [16] can be an advantage to study the scattering using an asymmetric box; in fact, this allows to access 
more non-degenerate low-momentum modes for a given volume. 

Because we are are interesting in resonant scattering, let us introduce a new particle in the system with mass 
related to the physical mass of n by 2m K < m§ < 4m K ; at the moment, it is not interacting with the particles n. In 
Figure 1 (Left) the energy level is the horizontal line; it intersects the two-particle levels at various system sizes L. 
If we introduce a three point interaction between <j) and n, in a Minkowski space and in an infinite volume, we can say 
the particle <p becames unstable and it decays into two ^-particles. In Euclidean space and in a finite volume we can 
only say that, because of the interaction, the energy eigenstates are a mixture of <j) and 2%. This is manifested as ALCs 
in the energy levels as shown in Figure 1 (Right). 

When the interaction is turned on Eq. (5) is still valid, but the momentum quantization condition is replaced with 
the Liischer formula which connects the momentum p = \p\ with the scattering phase shift dr. 



Sl{p) = mod ^; 1 = 



(6) 
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FIGURE 1. (Left) The spectrum of a system of two non-interacting particles of mass m n = 0.468; the horizontal line describes 
the particle <j> at rest of mass = 1.405. With these parameters the intersection between <j> and the two-particle level n 2 = 1 is set 
at L = 12. (Right) Avoided level crossings where on the Left there were intersections between (j> and In. 



where (q) is determined when I = by 

tan0(g) = - rc 3/ \ , 0(0) =0. (7) 

The function Zoo(s;<? 2 ) is the generalized zeta function defined in [3]. The determination of the resonance parameters 
follows immediately. Using Monte Carlo simulations we can determine the two-particle energy spectrum (a plot 
looking like Figure 1 (Right)). Then employing Eq. (5) for each value of E(L) we can determine the corresponding 
value of p(L) and from Eq. (6) the value of 8(E). Near the resonance we can fit 8(E) to the relativistic Breit-Wigner 
formula 




and therefore we can determine the two parameters m§ and r^. Unfortunately, the method described up to now gives 
poor results when the width of the resonance is large (see the A resonance case, analyzed in [17]). In this case the 
ALCs are washed out, i.e. we cannot locate the presence of the resonance in the way possible using data like that 
shown in Figure 1 (Right). 



Probability distribution method 

An alternative method to solve the previous problem is based on a different way to analyze the finite volume energy 
spectrum [15]. The basic idea is to construct the probability distribution W(p) according to the prescriptions: 

1 . Measure the two-particle spectrum E n (L) for different values of L and for n = 1 , • • • ,N; 

2. Determine p n (L) using Eq. (5), for each n and L; 

3. Choose a suitable momentum interval \pi,P2] and introduce an equal-size momentum bin with length Ap; 

4. Count how many eigenvalues p n (L) are contained in each bin; 

5. Normalize this distribution in the interval \p\,pi\. 

It is possible to show that the probability distribution W(p) is given by W(p) = cZn=\ \p'n( L )] and differentiating 
Eq. (6) with respect to L, it turns out: 



p„=i L <t>'(qn(p)) 



[c is a normalization constant]. (9) 



Let us introduce Wq(p), determined by Eq. (9) with 8(p) = and L n (p) corresponding to the free energy levels. 
The authors of [15] showed that in order to subtract the background (free n particles) it is convenient to consider the 



subtracted probability distribution W(p) =W(p)— Wq(p). In the infinite volume limit, this last quantity is determined 
by 8(p) alone and close to the resonance, assuming a smooth dependence on p for the other quantities, it follows a 
Breit-Wigner shape with the same width. The main task of our work is to test this method on an effective field theory 
where a resonance emerges and then exploit the qualitative and quantitative knowledge in the real case of QCD. 

A SIMPLE MODEL OF SCALAR MESON DECAY 

We examine a non-linear sigma model, characterized by a field E G SU(2), coupled to a scalar field by a three-point 
interaction. The action is given by 

S = j >x{/?Tr(d M Z^4) -ii£j3Tr + + ^<j) x (- + m 2 ^ + X(j) x p Tr (d^d^l + v)^ , (10) 

where the relation between the field E and the three (pion) fields n, is given by Y. x = exp ■ The coupling j3 is 

related to the pion decay constant / by j3 = f 2 /4. The coefficient v is introduced and tuned so that (<j>) = 0. The action 
is invariant under global chiral transformation E — > UrLU^, where Ul,Ur G SU(2), apart from an explicit symmetry 
breaking term m^j3Tr(E x + E^). It is instructive to consider the action for j3 — > °°, i.e. when the three pions do not 
interact with each other. In this limit, 

S = \ J d^x{Ki{-d^+ml)ni + ^ x {-d^+ml)^ x + X(^ x d^nid^ni + Xf 2 v^ x } , (11) 

and we can see clearly the characteristic coupling term where the derivatives of the two pions appear. When we 
discretize this theory on the lattice we note it is possible to update both fields <p and E by means of the heatbath 
algorithm [18]. Naturally, this model closely resembles a linear sigma model, although the (pnn interactions differ 
slightly. A "bottom up" construction was used to ensure similar prescriptions could be followed to couple mesons with 
general spin. An example of this for the p -meson is found in [5]. 

MONTE CARLO SIMULATION 

In order to determine the two-particle spectrum we first introduce the partial Fourier transform (PFT) of the field %: 

1 — 2n 

n(n,t) = - Yn{x,t)e ap , Pi = — ni , n, ■•= 0,- • • ,L, 1 . (12) 

V x Li 

Then we study operators with zero total momentum and zero isospin: 

3 

O n {t) = ^ i {n 1 t)% i {-n,t); (13) 

i=l 

in particular we take in account five different operators, corresponding to n 2 = 0, 1 , 2, 3 , 4. A sixth operator, that clearly 
has the correct quantum number, is the PFT of the field <p with p = 0. To determine the energy levels we use a method, 
introduced in [8], based on a generalized eigenvalue problem applied to the correlation matrix function Cy (?) = (OiOj), 
i.e. a matrix whose elements are all possible correlators between the six operators. An example of spectrum is given 
in Figure 2 where the dashed lines represent the free two-pion spectrum. In this plot the parameters have been chosen 
to show an ALC between the level and the 2% level for n 2 = 1. The gap between these two levels is small compared 
to the difference between the others (the width of the resonance is small) and therefore it is simple to determine the 
position of the resonance, i.e. in this case the ALC is not washed out. 

Unfortunately, it looks like it is not possible to tune the parameters of this model to increase the width of the 
resonance arbitrarily. Two different problems arise, depending on the value of j3. If j3 < 0.25, a saturation effect 
appears so even if we increase the value of X, the width of the resonance does not increase. On the other hand if j3 is 
big, a vacuum instability effect prevents simulations beyond a X max value. At the moment, we are trying to understand 
the origin of these problems in more detail through perturbative calculations. 




FIGURE 2. The two-pion spectrum. The ALC is seen around the n 2 
two-pion spectrum, take in account the discretization artefacts. 



1 two-pion state. The dashed lines, which show the free 



CONCLUSION AND OUTLOOK 

The problem of studying the width of a resonance on the lattice has attracted much attention recently. The main tool 
to get results is based on the study of the volume dependence of the two-particle states spectrum. A new method to 
analyze the data, based on the concept of probability distribution, has been proposed recently. In this work we have 
started to study a lattice field model where the new method will be tested. This model is important not only from a 
qualitative point of view, but it will be quantitatively important: we can tune the parameters of models like this to 
obtain correlation functions with the same precision of a QCD Monte Carlo study and therefore predict the precision 
with which we can obtain the QCD resonance parameters. Investigations of the scalar glueball following a prescription 
like this have begun. 
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